Light Quark Masses in Multi-Quark Interactions 
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We suggest and discuss in detail a multi-quark three flavor Lagrangian of the Nambu - Jona- 
Lasinio type, which includes a set of effective interactions proportional to the current quark masses. 
It is shown that within the dynamical chiral symmetry breaking regime, the masses of the pseudo 
Goldstone bosons and their chiral partners, members of the low lying scalar nonet, are in perfect 
agreement with current phenomenological expectations. The role of the new interactions is analyzed. 
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I. INTRODUCTION 

A long history of applying the Nambu - Jona-Lasinio 
(NJL) model in hadron physics shows the importance of 
the concept of effective multi-quark interactions for mod- 
elling QCD at low energies. Originally formulated only in 
terms of four-fermion chiral-symmetric couplings [l|, , 
the model has been extended to the realistic three flavor 
and color case with U(1)a breaking six-quark 't Hooft 
interactions [3T4l6j and a set of eight-quark interactions 
[ijj ■ The last ones complete the number of vertices which 
are important in four dimensions for dynamical chiral 
symmetry breaking [l8|> Gil • 

The explicit breaking of chiral symmetry in the model 
is described by th e q uark mass term of the QCD La- 
grangian, e.g. 0, [Hf . As a result, deviations from the 
exact symmetry predictions are expressed by functions of 
the light quark masses. The current quark mass depen- 
dence is of importance for several reasons, in particular 
for the phenomenological description of meson spectra 
and meson-meson interactions, and for the critical point 
search in hot and dense hadronic matter. In the latter 
case it has a strong impact on the phase diagram. The 
mass effects may lead to a different phase structure. For 
instance, the large mass difference between s and u(d) 
quarks may disfavor the formation of the color-flavor- 
locked phase at intermediate density, and the conjecture 
regarding the two critical points structure finally may not 



be true [22 1 



The explicit chiral symmetry breaking (ChSB) by the 
standard mass term of the free Lagrangian is only a part 
of the more complicated picture arising in effective mod- 
els beyond leading order [23] . Chiral perturbation theory 
[24l - |27| gives a well-known example of a self consistent ac- 
counting of the mass terms, order by order, in an expan- 
sion in the masses themselves. In fact, NJL-type models 
should not be an exception from this rule. If one con- 
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siders multi-quark effective vertices, to the extent that 
't Hooft and eight-quark terms are included in the La- 
grangian, certain mass dependent multi-quark interac- 
tions must be also taken into account. It is the purpose 
of this paper to study such higher order terms in the 
quark mass expansion. In particular, we show the ability 
of the model with new quark-mass-dependent interac- 
tions to describe the spectrum of the pseudo Goldstone 
bosons, including the fine tuning of the rj-r]' splitting, and 
the spectrum of the light scalar mesons: <t(600), k(850), 
/ (980), and a (980). 

There are several motivations for this work. In the 
first place, the quark masses are the only parameters of 
the QCD Lagrangian which are responsible for an ex- 
plicit ChSB, and it is important for the effective theory 
to trace this dependence in full detail. In this paper it 
will be argued that it is from the point of view of the 
1/N C expansion that the new quark mass dependent in- 
teractions must be included in the NJL-type Lagrangian 
already when the U(1)a breaking 't Hooft determinantal 
interaction is considered. This point is somehow com- 
pletely ignored in the current literature. 

A second reason is that nowadays it is getting clear 
that the eight-quark interactions, which are almost 
inessential for the mesonic spectra in the vacuum, can 
be important for the quark matter in a strong magnetic 
background [2S--32] . We will show that there is a set 
of the effective quark-mass-dependent interactions which 
are of importance here and have not been considered yet. 

A further motivation comes from the hadronic matter 
studies in a hot and dense environment. It is known that 
lattice QCD at finite density suffers from the numerical 
sign problem. Thus, the phase diagram is notoriously 
difficult to compute "ab initio" , except for an extremely 
high density regime where perturbative QCD methods 
are applicable. In such circumstances effective models 
designed to shed light on the phase structure of QCD 
are valuable, especially if such models are known to be 
successful in the description of the hadronic matter at 
zero temperature and density. Reasonable modifications 
of the NJL model are of special interest in this context 
and our work aims at future applications in that area. 
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II. EFFECTIVE MULTI-QUARK 
INTERACTIONS 



The chiral quark Lagrangian has predictive power for 
the energy range which is of order A ~ 4^/^ ~ 1 GeV 
[33l j . A characterizes the spontaneous chiral symmetry 
breaking scale. Consequently, the effective multi-quark 
interactions, responsible for this dynamical effect, arc 
suppressed by A, which provides a natural expansion 
parameter in the chiral effective Lagrangian. The scale 
above which these interactions disappear and QCD be- 
comes perturbative enters the NJL model as an ultravi- 
olet cut-off for the quark loops. Thus, to build the NJL 
type Lagrangian we have only three elements: the quark 
fields q, the scale A, and the external sources x^ which 
generate explicit symmetry breaking effects - resulting in 
mass terms and mass-dependent interactions. 

The color quark fields possess definite transformation 
properties with respect to the chiral flavor U(3)l x U(3)r 
global symmetry of the QCD Lagrangian with three 
massless quarks (in the large N c limit). It is con- 
venient to introduce the U(3) Lie-algebra valued field 
£ = (s a - ip a )\K, where s a = qX a q, p a = qX a ij 5 q, 
and a = 0, 1, . . . , 8, Ao = \/2/3 x 1, A a being the stan- 
dard ,577(3) Gell-Mann matrices for 1 < a < 8. Un- 
der chiral transformations: q' — VbMb. + Vl<7l, where 
qn = Pnq,qL = Phq, and P R>L = |(1 ± 75). Hence, 
£' = VrEVI, and £+' = VlE+V^. The transformation 
property of the source is supposed to be x' — VrxVl ■ 

Any term of the effective multi-quark Lagrangian with- 
out derivatives can be written as a certain combination 
of fields which is invariant under chiral SU(3)r x SU (3) l 
transformations and conserves C, P and T discrete sym- 
metries. These terms have the general form 



(1) 



where g~i are dimensionless coupling constants (starting 
from eq. (f2~Tj) the dimensional couplings gi — <?i/A 7 will 
be also considered). Using dimensional arguments we 
find a + 3/3 — 7 = 4, with integer values for a, /3 and 7. 

We obtain a second restriction by considering only the 
vertices which make essential contributions to the gap 
equations in the regime of dynamical chiral symmetry 
breaking, i.e. we collect only the terms whose contribu- 
tions to the effective potential survive at A -> 00. We get 
this information by contracting quark lines in Li , finding 
that this term contributes to the power counting of A in 
the effective potential as ~ A 2 ^ -7 , i.e. we obtain that 
2/3 — 7 > (we used the fact that in four dimensions 
each quark loop contributes as A 2 ). 

Combining both restrictions we come to the conclusion 
that only vertices with 



a + (3 < 4 



(2) 



must be taken into account in the approximation consid- 
ered. On the basis of this inequality one can conclude 



that (i) there are only four classes of vertices which con- 
tribute at a — 0; those are four, six and eight-quark in- 
teractions, corresponding to (3 — 2, 3 and 4 respectively; 
the (8 = 1 class is forbidden by chiral symmetry require- 
ments; (ii) there are only six classes of vertices depend- 
ing on external sources Xj they are: a = 1,0 = 1,2,3; 
a = 2, (3 = 1, 2; and a = 3, /3 = 1. 

Let us consider now the structure of multi-quark ver- 
tices in detail. The Lagrangian corresponding to the case 
(i) is well known 



L mt = ^tr(S t E) + ^(detE + detS t ) 
+ % (trSt S ) 2 + ||tr(l]t SS t s) . 

It contains four dimensionful couplings G, k, gi, gi- 
The second group (ii) contains eleven terms 



(3) 



where 



in 

i=0 



L = -tr(£t x + x t E ) 

L\ = —-^eijuemnl^imXjnXkl+h.C. 

7/2 — 'J^^ijk^rnnlXiyn^-'jrL^-'kl ~t~ h.C. 

L 3 = |§tr (£tEEt x ) + h.c. 

L 4 = Htr(EtE)tr(Et x )+fc.c. 

L 5 = ||tr (St x E+ X ) + h.c. 

L 6 = gtr(££t xx t + E t Ex t x ) 

L 7 = VL(tr^ X + h.c.) 2 



(4) 



7* 



7, 



A4 (tr£t x -/i.c.)' 

"tr(x t x)tr(x t S)+ft.c (5) 



7 9 = --r|tr (E + xx f x) +h.c 



A 2 
A 2 



Each term in the Lagrangian 7g is hermitian by itself, but 
because of the parity symmetry of strong interactions, 
which transforms one of them into the other, they have 
a common coupling gg. 

Some useful insight into the Lagrangian above can be 
obtained by considering it from the point of view of the 
1 /N c expansion. Indeed, the number of color components 
of the quark field q l is N c , hence summing over color 
indices in E gives a factor of N c , i.e. one counts E ~ N c . 

The cut-off A that gives the right dimensionality to the 
multi-quark vertices scales as A ~ N® — 1 . On the other 
hand, since the leading quark contribution to the vacuum 
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energy is known to be of order N c , the first term in ([3j) 
is estimated as N c , and we conclude that G ~ 1/N C . 

Furthermore, the U(1)a anomaly contribution (the 
second term in ([3])) is suppressed by one power of 1/N C , 
it yields k ~ l/iV*. 

The last two terms in ([3]) have the same N c count- 
ing as the 't Hooft term. They are of order 1. Indeed, 
Zweig's rule violating effects are always of order 1/N C 
with respect to the leading order contribution ~ N c . 
This reasoning helps us to find g\ ~ l/-^c- The term 
with g 2 ~ is also 1/N C suppressed. It represents 

the next to the leading order contribution with one in- 
ternal quark loop in N c counting. Such vertex contains 
the admixture of the four-quark component qqqq to the 
leading quark- antiquark structure at N c — > oo. 

Next, all terms in eq. ([3]), except Lq, are of order 1. 
The argument is just the same as before: this part of 
the Lagrangian is obtained by succesive insertions of the 
X- held (x counts as x ~ 1) m place of S fields in the 
already known 1/N C suppressed vertices. It means that 
«i, 59,ffio ~ ^2, 35, 36, 37, 3s ~ 1/^V C 2 , and g 3 ,g 4 ~ 

There are two important conclusions here. The first 
is that at leading order in 1/N C only two terms con- 
tribute: the first term of eq. ([3]), and the first term 
of eq. ([5]). This corresponds exactly to the standard 
NJL model picture, where mesons are pure qq states. 
At the next to leading order we have thirteen terms ad- 
ditionally. They trace the Zweig's rule violating effects 
(k, «i, K2, 51,54,57, 38, 5io), and an admixture of the four- 
quark component to the qq one (32,33,35,36, 39)- Only 
the phenomenology of the last three terms from eq. ([3J 
has been studied until now. We must still understand 
the role of the other ten terms to be consistent with the 
generic 1/N C expansion of QCD. 

The second conclusion is that the N c counting justifies 
the classification of the vertices made above on the basis 
of the inequality @. This is seen as follows: the equiva- 
lent inequality \(a + (3)/2~\ < 2 is obtained by restricting 
the multi-quark Lagrangian to terms that do not vanish 
at N c -> 00 (it follows from dJ) that (3 - [7/2] > by 

noting that gi ~ l/ivj 7 ^, where [7/2] is the nearest 
integer greater than or equal to 7/2). 
The total Lagrangian is the sum 

L = qiy^d^q + L int + L x . (6) 

In this SU(3)l x SU(3)r symmetric chiral Lagrangian 
we neglect terms with derivatives in the multi-quark in- 
teractions, as usually assumed in the NJL model. We 
follow this approximation, because the specific questions 
for which these terms might be important, e.g. the radial 
meson excitations, or the existence of some inhomoge- 
neous phases, characterized by a spatially varying order 
parameter, are not the goal of this work. 

Finally, having all the building blocks conform with 
the symmetry content of the model, one is now free to 
choose the external source x- Putting x = -M/2, where 



we obtain a consistent set of explicitly breaking chiral 
symmetry terms. This leads to the following mass de- 
pendent part of the NJL Lagrangian 

8 

L x -> Lfj, = -qmq + L ■ (7) 

where the current quark mass matrix m is equal to 

m = M + ^-(detM)M- 1 + ^M 3 
A v ' 4A 2 

+ g| (trM 2 ) ^, (8) 

and 

^-"1 = 2~K3^ijkCrnnl-Mi m Y!ij n Yifci + h.C. 

L' 3 = ^tr (StEEtM) + h.c. 

L' 4 = J^tr (StE) tr (tfM) + h.c. 

L' 5 = ^tr (tfMTtM) + h.c. (9) 

L^ = ^tr [M 2 (EE* + E^E)] 

L' 7 = ^(tr^M + h.c.f 

L> s = $-(trtfM-h.c.) 2 

Let us note that there is a definite freedom in the def- 
inition of the external source \. In fact, the sources 

X (Cl) = X + j- (detx^xixh)' 1 + -^XX f X 

+ §tr(x f x)x (10) 

with three independent constants c$ have the same sym- 
metry transformation property as x- Therefore, we could 
have used x^ Ci ' everywhere that we used x- As a result, 
we would come to the same Lagrangian with the follow- 
ing redefinitions of couplings 

_/ _ , Cl _ _, _ _ 
Kl Kl = Kl + _ g 5 _> = g 5 - K 2Clj 

-I - , «2 - -/ - , K 2 

37 -> 3 7 = 37 + yci, 3s -> g 8 = g s + —c x , 
39 -> g'g = 39 + c 2 - 2k\c x , 

5io -> 5io = 5io + C3 +2rtiCx. (II) 

Since Ci are arbitrary parameters, this corresponds to a 
continuous family of equivalent Lagrangians. This family 
reflects the known Kaplan - Manohar ambiguity [3J-[33| 
in the definition of the quark mass, and means that sev- 
eral different parameter sets (fTTj) may be used to repre- 
sent the data. In particular, without loss of generality 
we can use the reparametrization freedom to obtain the 
set with = g'g = g[o = 0. 
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The effective mufti-quark Lagrangian can be written 
now as 



L = qiWdp - m)q + L mt + V L \ 



i=2 



(12) 



It contains eighteen parameters: the scafe A, three pa- 
rameters which are responsibie for explicit chiral sym- 
metry breaking (i u , f*d, fts, and fourteen interaction cou- 
plings G,k,Ri,R 2 , gi,...,g w . Three of them, Ki,g s ,g w , 
contribute to the current quark masses m. Seven more 
describe the strength of multi-quark interactions with ex- 
plicit symmetry breaking effects. These vertices contain 
new details of the quark dynamics which have not been 
studied yet in any NJL-type models. 



III. FROM QUARKS TO MESONS: 
STATIONARY PHASE CALCULATIONS 



The model can be solved by path integral bosonization 
of this quark Lagrangian. Indeed, following [7J we may 
equivalently introduce auxiliary fields s a = qX a q, p a = 
<?*75A a 9, and physical scalar and pseudoscalar fields a = 
CaA a , </> — 4>a A a • In these variables the Lagrangian is a 
bilinear form in quark fields (once the replacement has 
been done the quarks can be integrated out giving us the 
kinetic terms for the physical fields <f> and a) 

L = q [i^dn - (a + 175^)] q + L auxi 
Laux = s a a a + p a 4> a - s a m a + L int (s,p) 



8 
i=2 



(13) 



It is clear, that after the elimination of the fields a, 4> 
by means of their classical equations of motion, one can 
rewrite this Lagrangian in its original form (1121) . On the 
other hand, written in terms of auxiliary bosonic vari- 
ables, the Lagrangian becomes 

Li nt (s,p) = L Aq + L 6q + L^j + L^j, 



G 

2A 2 " 

R 

4A 



Liq(s,p) = ^ (sI+pI) , 



Ls q (s,p) = -J^A a bcSa{sbSc - SpbPc 



(14) 



L^(s,p) = -^[dabed c de(s a Sb+PaPb){ScSd+PcPd) 
+ 4:fabefcdeS a S c p b Pd\ , 

and the quark mass dependent part is as follows 

= ■^^ A abc^a(SbSc-PbPc) , 



rl 94 I 2 1 2\ 

4 = 2A6 MbSfc ' Sa +Pa > ' 

L' 5 = -^J^^b^d (dabedede ~ fabefede) {s a S c ~ PaPc) , 

L' 6 = ■^j^HafJ.bdabedcdeiScSd+PcPd) , 

T I _ 97 / ^2 

1 ~ ^ aSa > ' 

T I 98 , -.2 

L 8 = -JZ {fJ-aPa) , 



(15) 



where 



-A-abc — ^&ijk&rri7il 

(A a )im(A&)jn(A c )fc;, (16) 



and the U (3) antisymmetric f a bc and symmetric d a bc con- 
stants are standard. 

Our final goal is to clarify the role of the mass- 
dependent terms described by the Lagrangian densites 
of eq. (fT5")) . We can gain some understanding of this by 
considering the low-energy meson dynamics which follows 
from our Lagrangian. For that we must exclude quark 
degrees of freedom in (|13[) . e.g., by integrating them out 
from the corresponding generating functional. The stan- 
dard Gaussian path integral leads us to the fermion de- 
terminant ! _which we expand by using a heat-kernel tech- 
nique [38l44l| . The remaining part of the Lagrangian, 
Laux, depends on auxiliary fields which do not have ki- 
netic terms. The equations of motion of such a static 
system are the extremum conditions 



dL -o dL -o 

ds a ' dp a 



(17) 



which must be fulfilled in the neighbourhood of the uni- 
form vacuum state of the theory. To take this into ac- 
count one should shift the scalar field a — > cr + M, 
The new cr-field has a vanishing vacuum expectation 
value (cr) = 0, describing small amplitude fluctuations 
about the vacuum, with M being the mass of constituent 
quarks. We seek solutions of eq. (|17l) in the form: 



l abc 



^abcMc 



Pa = ^ 2 6 WO^c + 



(18) 



Eqs. (IT71) determine all coefficients of this expansion giv- 
ing rise to a system of cubic equations to obtain h a , and 
the full set of recurrence relations to find higher order 
coefficients in (|18l) . We can gain some insight into the 
physical meaning of these parameters if we calculate the 
Lagrangian density L aU x on the stationary trajectory. In 
fact, using the recurrence relations, we are led to the re- 
sult 



La 



1 



1, 

2 ab 



(2). 



(19) 



— cr„ 
3 



h abc (T b(Tc 



(2) 



,(3) 



L' 



93 
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Ha[dabed c deSb{ScSd+ PcPd) - ZfabefcdePbPcSd] , sates, h9j} , h{\ 



From this one can see that h a define the quark conden- 



h\ b ' contribute to the masses of scalar and 
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pseudoscalar states, and higher order h's are the cou- 
plings that measure the strength of the meson-meson in- 
teractions. 

We proceed now to explain the details of determining 
h. In the following only the first coefficients h a , h^ b , and 

^ab wm ^ e °^ m t eres t to us. In particular, eq. (fTT)) states 
that h a = 0, if a ^ 0, 3, 8, while h a (a = 0, 3, 8), after the 
convenient redefinition to the flavor indices i = u,d,s 



kr, = 



1 



2V3 



y/2 V2 y/2 ' 
n/3 -V3 



(20) 



I 1 -2 
satisfy the following system of cubic equations 



A, + ^Ujkhjhk + ^ (2G + 5 i h 2 + g Af ih) + ^-hf 



-j [333/if + 34 h 2 + 2(.g 5 + gojmhi + 437 -fih] 

+K2tijkHjhk = 0. 



(21) 



Here A* = Mi — mi ; t ijk is a totally symmetric quantity, 
whose nonzero components are t ue i s = 1 ; there is no sum- 
mation over the open index i but we sum over the dummy 
indices, e.g. h 2 = h 2 t + h 2 d + h 2 , jj,h = fi u h u + ^ d h d + n s h s . 
In particular, eq. (|8|) reads in this basis 



M , S9 2 , JlO 2 

mi = ^ ( 1 + — ^ + —11 



ijk ft j ftk ■ 



(22) 



For the set 39 = 310 = k± = the current quark mass rrii 
coincides precisely with the explicit symmetry breaking 
parameter /x,. 

Note that the factor multiplying hi in the third term 
of eq. (f2"Tj) is the same for each flavor. This quantity also 
appears in all meson mass expressions, and there is no 
further dependence on the couplings G, 31 , 34 involved for 
meson states with a = 1, . . . , 7. Thus there is a freedom 
of choice which allows to vary these couplings, conden- 
sates and quark masses /i,-, without altering this part of 
the meson mass spectrum. 

It is now straightforward to obtain the inverse matrices 
to hty and , namely 

-2 ( h ab) 1 = (2G + 31 h 2 + 34/z/i) 5 ab + AgxKhi 
+3A abc [nh c + 2n 2 p c ) + g2h r h c (d a bed C re + 2d ace d bre ) 

^dSf^r^c (^afte^cre ~t~ d ace d bre -{- d are d bce ^ 

+ 234 (f+ahb + fJ-bha) + g^rV-c (daredbce ~ farefbce) 

+g6HrfJ-cdabed C re + ^97^a^b- (23) 

-2 (h$) 1 = (2G + 9l h 2 + g^h) S ab 
-3A abc (nh c + 2k 2 Hc) + 92h r h c (d abe d cre + 

^farefbce^) 

J rgzft r h c (dated 

ere ~l~ farefbce facefbre) 

-g 5 fj, r p c (d 

are dbce farefbce ) 
+geHr^cdabed C re ~ 4:g S fl a fl b . (24) 

These coefficients are totally defined in terms of h a and 
the parameters of the model. 



IV. FROM QUARKS TO MESONS: HEAT 
KERNEL CALCULATIONS 



We now turn our attention to the total Lagrangian of 
the bosonized theory. To write down this Lagrangian we 
should add the terms coming from integrating out the 
quark degrees of freedom in (TIB")) to our result (TIT)]) . For- 
tunately, the result is known. One can find all necessary 
details of such calculations for instance in 
used the modified heat kernel technique 



38 , where we 
MEI devel- 



oped for the case of explicit chiral symmetry breaking. 
Here we quote the main outcome. The a tadpole term 
must be excluded from the total Lagrangian. This gives 
us a system of gap equations. 



hi + ^Mi [3/ 



(3M 2 - M 2 ) h] = 0. 



(25) 



Here N c — 3 is the number of colors, and M 2 = M 2 + 
M 2 + M 2 . The factors h (i = 0, 1, . . .) are the arith- 
metic average values h = |[Jj(M 2 ) + Ji{Mj) + J^M 2 )], 
constructed from the one-quark-loop integrals 



Ji(m 2 ) 



dt 
T 2 ^ 



p{tk 2 )e 



(26) 



with the Pauli-Villars regularization kernel [42J, |43| 

P {tk 2 ) = 1 - (1 + iA 2 ) exp(-iA 2 ). (27) 
In the following we need only to know two of them 



A 2 

J (w 2 ) = - m 2 In ( 1 + — 

mr 



and 



A 2 

Ji(m?) = In I 1 + 

mr 



A 2 



A 2 



(28) 



(29) 



From now on we will consider the case with an exact 
SU(2) isospin symmetry, i.e. jj, u = pd = A 7^ ^s, and 
M u = Md = M ^ M s . We also restrict ourselves to 
small perturbations, so we retain terms in the bosonized 
Lagrangian which are quadratic in the perturbations <j) 
and a. To this order we obtain 



NJi 
16tt 2 
N c h 



tr [(d^a) 2 + (O^) 2 



NJ 



4tt 2 



2(2M 2 + M 2 )a 2 + (M 2 + 5M 2 )a 2 



+ {7M 2 - M 2 )af + (M + M S )(M + 2M s )a) 



(M s - M)(2M S ~ M)q> 2 } 



lib 



V a Ob 



I Kb 



(2). 



h + 



(30) 
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where A n? = M 2 - M 2 , # = £^0?, $ = £/=4$- 
The kinetic term requires a redefinition of meson fields, 

47T 2 

""a = Serf, a = ff 2 = -rf-f-, (31) 

to obtain the standard factor 1/4. The Lagrangian (|30|) 
in the chiral limit, m — 0, leads to the conserved vector, 
and axial- vector, currents. The matrix elements 
of axial- vector currents 

(0|^(0)|^(p)> = ip^/ 06 (32) 

define the weak and electromagnetic decay constants of 
physical pseudoscalar states (see details in 38]). In 
fact, we obtain that all new information about the mass- 
dependent interactions is explicitly absorbed in the last 

(1 2) 

two terms of the Lagrangian, where the matrices h i 
are block diagonal and mix only in the (0,8) sector, see 
eqs. ([23]) and (j24|) . There is also an implicit dependence 
through the gap and stationary phase equations. 

V. FIXING PARAMETERS 

Now let us fix the values of the various quantities in- 
troduced. After choosing the set k,\ = gg = g w = we 
still have to fix fourteen parameters: h.,rh,m s ,G, k, k% 
and gi , . . . , <78 ■ Note that there are two intrinsic restric- 
tions of the model, namely, the stationary phase (|21|) and 
the gap (|25|) equations, which must be solved selfconsis- 
tently. This is how the explicit symmetry breaking is 
intertwined with the dynamical symmetry breaking and 
vise versa. We use ([25]) to determine h, h s through A, M s 
and M. The ratio M s /M is related to the ratio of the 
weak decay constants of the pion, /„■ = 92 MeV, and the 
kaon, Jk = 113 MeV. Here we obtain 

^ = 2^- 1 = 1.46. (33) 
M U 

Furthermore, the two eqs. (|2"TT) can be used to find the 
values of A and M if the parameters to, to s , G, k, K2, 
gi,...,gt are known. Thus, together with gg, we have 
at this stage thirteen couplings to be fixed. Let us con- 
sider the current quark masses to and m s to be an input. 
Their values are known, from various analyses of the chi- 
ral treatment of the light pseudoscalars, to be around 
in = 4 MeV and m s = 100 MeV 0]. Then the remain- 
ing eleven couplings can be found by comparing with 
empirical data. One should stress the possibility (which 
did not exist before the inclusion of mass-dependent in- 
teractions) to fit the low lying pseudoscalar spectrum, 
niK = 138 MeV, m K = 494 MeV, m v = 547 MeV, 
m^i = 958 MeV, the weak pion and kaon decay constants, 
f„ = 92 MeV, f K = 113 MeV, and the singlet-octet mix- 
ing angle 9 P = — 15° to perfect accuracy. One can deduce 
that the couplings k 2 and gg are essential to improve the 
description in the pseudoscalar sector; in particular, gg 
is responsible for fine tuning the Tj—rf mass splitting. 



The remaining five conditions are taken from the scalar 
sector of the model. Unfortunately, the scalar channel in 
the region about 1 GeV became a long-standing problem 
of QCD. The abundance of meson resonances with ++ 
quantum numbers shows that one can expect the pres- 
ence of non-qq scalar objects, like glueballs, hybrids, mul- 
tiquark states and so forth [45| ■ This creates known diffi- 
culties in the interpretation and classification of scalars. 
For instance, the numerical attempts to organize the t/(3) 
quark-antiquark nonet based on the light scalar mesons, 
a or / Q (600), a (980), k(850), / (980), in the framework 
of NJL-type models have failed (see, e.g. [IM S 0- 
l48l|). The reason is the ordering of the calculated spec- 
trum which typically is m a < m ao < m K < to/ , as 
opposed to the empirical evidence: m K < m aa ~ to/ . 

On the other hand, it is known that a unitarized non- 
rclativistic meson model can successfully describe the 
light scalar meson nonet as qq states with a meson-meson 
admixture [49(. Another model which assumes the mix- 
ing of qq-states with others, consisting of two quarks and 
two antiquarks, q 2 q 2 (Hoj . yields a p ossible description of 
the ++ meson spectra as well 15 1[ 15211 . The well known 
model of Close and Tornqvist [53| is also designed to 
describe two scalar nonets (above and below 1 GeV). 
The light scalar nonet below 1 GeV has a core made 
of q 2 q 2 states with a small admixture of a qq compo- 
nent, rearranged asymptotically as meson-meson states. 
These successful solutions seemingly indicate on the im- 
portance of certain admixtures for the correct description 
of the light scalars. Our model contains such admixtures 
in the form of the appropriate effective multi-quark ver- 
tices with the asymptotic meson states described by the 
bosonized qq fields. We have found, that the quark mass 
dependent interactions can solve the problem of the light 
scalar spectrum and these masses can be understood in 
terms of spontaneous and explicit chiral symmetry break- 
ing only. Indeed, one can easily fit the data: to ct = 600 
MeV, m ao = 980 MeV, m K = 850 MeV, m fo = 980 MeV 
with the input value g 2 = 0. In this case we obtain for 
the singlet-octet mixing angle 8 S roughly 9 S = 19°. 

To many readers our success with scalars may seem 
trivial: we have five parameters to fit five numbers. What 
is not trivial, however, is that the overall result of the fit 
is also in an agreement with phenomenological expec- 
tations. To compare, if we try instead to fit the sec- 
ond scalar nonet f (1370), a (1450), ^(1430), f (1500) 
with the same input, our attempt fails. The best that 
we can do is the values to/ = 1220 MeV, m ao = 1406 
MeV, m K * = 1506 MeV, m/< = 1786 MeV. However, 
even these unreasonable masses come out only together 
with the very large ratio m s /rh = 36 and phenomcnolog- 
ically unacceptable values for constituent quark masses 
M = 631 MeV and M s = 919 MeV. 

We obtain and understand the empirical mass assign- 
ment inside the light scalar nonet as a consequence of 
the quark-mass dependent interactions, i.e. as the re- 
sult of some predominance of the explicit chiral symme- 
try breaking terms over the dynamical chiral symmetry 
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TABLE I: Parameter sets of the model: m, m s , and A are given in MeV. The couplings have the following units: [G] = GeV 2 , 
[k] = GeV -5 , [gi] = [32] = GeV -8 . We also show here the values of constituent quark masses M and M B in MeV. 



Sets 


m 


m s 


M 


M s 


A 


G 


—k 


9i 


92 


a 


4.0* 


100* 


361 


526 


837 


8.96 


93.0 


1534 


0* 


b 


4.0* 


100* 


361 


526 


837 


7.06 


93.3 


3420 


0* 



TABLE II: Explicit symmetry breaking interaction couplings. The couplings have the following units: [m] = GeV 1 , [^2] — 
GeV -3 , [g 3 ] = [ 54 ] = GeV -6 , [g 5 ] = [g e ] = [57] = [gs] = GeV -4 , [g 9 ] = [g 10 ] = GeV -2 . 



Sets 




K2 


-53 


-94, 


95 


-9o 


-97 


58 


99 


9io 


a 


0* 


9.05 


4967 


661 


192.2 


1236 


293 


52.2 


0* 


0* 


b 


0* 


9.01 


4990 


653 


192.5 


1242 


293 


51.3 


0* 


0* 



breaking ones for these states. Indeed, let us consider the 
difference 

m 2 a — m 2 = 2q 2 ( — -) 

y \H ao Hj 

- 2 (M s + 2M)(M S - M). (34) 

The sign of this expression is a result of the competition 
of two terms. In the chiral limit both of them are zero, 
since at /*, ji s = we obtain M — M s and H ao = H K , 
for H ao and H K being positive. The splitting H K > H aQ 
is a necessary condition to get m ao > m K . The following 
terms contribute to the difference 

H K - H ao = K(h s - h) + 2k 2 (Ms - A) - 92{h 2 s + hh s - 2h 2 ) 
+ y \ 2nsh s + Hs 'h + jj,h s - iftfij 

+ gsKp-s - A) + -| {A - A 2 ) ■ (35) 

Accordingly, from this formula we deduce the "anatomy" 
of the successful numerical fit: 

"4 - ml = ([0.007]„ + [0.076] K2 + [0] ff2 

+ [0.832] 93 + [0.003] S5 + [-0.269] 96 

- [0.41] M = 0.24) GeV 2 , (36) 

where the contributions of terms with corresponding cou- 
pling (see eq. ([33]) ) are indicated in square brackets. The 
last number, marked by M, is the value of the last term 
from It is a contribution due to the dynamical chiral 
symmetry breaking (in the presence of an explicit chiral 
symmetry breaking) . One can see that the (73-interaction 
is the main reason for the reverse ordering m ao > m R , 
the coupling g§ being responsible for the fine tuning of 
the result. 

Let us now show the result of our global fitting of the 
model parameters. We collect them in two tables. Two 
sets (a) and (b) are shown. The difference is the fit- 
ted value of the a mass: in (a) m a — 600 MeV, in (b) 
m a = 500 MeV. Table 1 contains the standard set of pa- 
rameters, which are known from previous considerations. 
Their values are not much affected by the quark mass 



effects. Table 2 contains the couplings which are respon- 
sible for the explicit chiral symmetry breaking effects in 
the interactions. Note that these couplings almost do not 
change from set (a) to (b). We have already learned (as 
seen again in Table 1) that higher values of g\ lead to the 
lower a mass (38[. This eight-quark interaction violates 
Zweig's rule, since it involves qq annihilation. The mixing 
angle 9 S is stable with respect to such changes, we obtain 
6 S = 19.4° in case (a), and S = 18.9° in case (b). The 
calculated values of quark condensates are the same for 
both sets: ~(uu)^ = 232 MeV, and -(ss)i = 206 MeV. 
Our calculated values for constituent quark masses agree 
with the ones found in 111 showing their insensi- 
tivity to the new mass-dependent corrections. 



VI. CONCLUDING REMARKS 

The purpose of this paper has been to take into ac- 
count the quark masses at next to leading order in the 
expansion of the effective multi-quark Lagrangian of the 
NJL-type. As a result a picture with some attractive 
new features has emerged. Let us summarize the details 
of such a picture. 

The main qualitative difference between our result and 
previous calculations is the possibility to fit the low ly- 
ing pseudoscalar spectrum (the pseudo Goldstone h 
nonet) and weak decay constants of the pion and the 
kaon to perfect accuracy. The fitting of the 77 — 17' mass 
splitting together with the overall successful description 
of the whole set of low-energy characteristics is actually a 
solution for a long standing problem of NJL-type models. 
We expect that with such modifications the model is get- 
ting more appropriate not only for studying low-energy 
meson physics, but also in studies of the ground state of 
hadronic matter in an environment, which is known to 
be very sensitive to quark mass effects. 

With a set of new quark-mass dependent interactions 
we are also capable to describe the spectrum of the light 
scalar nonet. From that one can conclude that both spec- 
tra can be understood on the basis of the dynamical and 
explicit chiral symmetry breaking only. The splitting in- 
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side the scalar nonet is determined by two competing 
contributions: first it is due to the explicit symmetry 
breaking (embodied in the stationary phase part of the 
bosonized Lagrangian), second it is due to the dynami- 
cal symmetry breaking (see the heat kernel part of the 
bosonized Lagrangian). It is the first type of contribution 
that changes the ordering inside the light scalar nonet, 
as compared to the standard approach. 

Our result for the scalar sector, being promising by it- 
self, must be considered with some reservation. To report 
about a real success here, one should explain not only the 
mass spectrum of scalars, particularly the mass degener- 
acy of the /o(980) and ao(980) states (as we have done 



here), but answer some known challenges related with 
radiative decays of these states. Work in this direction is 
in progress. 
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